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Abstract 

The phenomenology of heavy quarkonia production in fixed target experiments and at the 
Tevatron and HERA colliders is reviewed. The latest theoretical results are presented and 
compared with data, with emphasis on the predictions of the factorization approach by Bodwin, 
Braaten and Lepage. 


* To appear, in a slightly shortened form, in the Proceedings of the XXXII Rencontres de Moriond, 
“QCD and High Energy Hadronic Interactions”, Les Arcs, March 1997. 




1 Introduction 


The production of heavy quarkonia has been subjected to intense study in the last two or 
three years, with tens of papers having being produced on the problem of J/0’s, y’s and T’s 
production in e + e“, yp, pp, pN, tiN collisions and also B decays. 

The reason for such a surge in interest was the appearance of a theoretical framework, the 
so called Factorization Approach (FA) by Bodwin, Braaten and Lepage |l|], which seems able 
to solve the theoretical problems that quarkonia production models faced in the past, and also 
to reconcile theoretical predictions with experimental data, previously in disagreement up to 
factors of fifty in some instances. 

In this talk I shall not review the Factorization approach in detail, leaving this theoretical 
introduction to other sources (see for instance b s a)- 1 shall also not discuss in detail the 
Color Singlet Model (CSM) || and the Color Evaporation Model (CEM)|| (the latter has also 
been recently compared to data and found able to describe at least some of them [j7. [|). 

I shall just recall how the factorization approach writes the quarkonium state H production 
cross section in the following form: 

a{ij -> QQ -*• H) -> QQ[n])(0 H (n)). (1) 

n 

According to this equation, the cross section for producing the observable quarkonium state H 
is factorized into two steps. In the short distance part a QQ pair of heavy quarks is produced in 
the spin/color/angular momentum state 2S+1 L^ = n by the scattering of the two light partons 
i and j. Successively this pair hadronizes into the quarkonium //, and (O h ( n)) is formally a 
Non Relativistic QCD (NRQCD) matrix element describing this non perturbative transition. 

An important feature of this equation is that also QQ pairs in a color octet state are allowed 
to contribute to the production of a color singlet quarkonium H\ their color is neutralized 
via a non perturbative emission of soft gluons. While the corresponding matrix elements are 
suppressed by the need of such an emission, the short distance coefficients can on the other 
hand be large, perhaps overcompensating the suppression of the non perturbative term. This 
explains why color octet contributions can play a very important role in predicting the total 
size of quarkonia production cross sections. 

Many theoretical items would be worth discussing about the Factorization Approach. Ac¬ 
tually, they would probably be worth a seminar (or more) by themselves. As I said, I shall 
however skip such a detailed discussion and rather concentrate on some selected phenomeno¬ 
logical outcomes of the theoretical investigations which have been carried on so far. I shall 
restrict myself to analyses of experimental data coming from pp collisions at the Tevatron, yp 
collisions at HERA and fixed target experiment, reviewing the results of these investigations. 
References will be provided to the theoretical papers, leaving to them the task of properly citing 
the experimental ones. 

The aim of the game will be to check whether the non perturbative matrix elements which 
can be extracted from these experimental data are mutually consistent with each other. In 
other words, we shall check whether the factorization approach can properly describe all the 
data with matrix elements truly universal and independent from the underlying short distance 
process, as they should. 
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2 The Tevatron data in pp collisions 


Beginning with the Tevatron data looks appropriate as the explanation of its anomalously 
large ip' production rate (a factor of fifty above the Color Singlet Model prediction) was the 
first phenomenological breakthrough of the FA. Braaten and Fleming [[JJ explained this large 
rate by assuming it was due to a color octet QQ pair originating via perturbative splitting 
from a large px gluon. While the non perturbative matrix element for such a color octet pair 
to produce a ip' is predicted by NRQCD to be about two order of magnitudes suppressed with 
respect to the one for a color singlet, the production rate of gluons (and hence of such pairs) 
is however so large that it can more than compensate for the suppression. Indeed, Braaten 
and Fleming could successfully describe the data fitting a value for the matrix element in good 
agreement with the theoretically predicted two-orders-of-magnitude suppression. 

This apparent success of the factorization approach on the ip' anomaly problem made imme¬ 
diately clear the potential importance of color octet mediated channels and stimulated similar 
research in other reactions: by the time of this Conference, the Braaten and Fleming’s paper 
has received more than 120 citations. 


Cross sections for large px J/ip’s and y’s production at the Tevatron have been analyzed, and 
a more detailed study of ip' has also been performed. It was found the theoretical curves could 
describe the shape of the data pretty well. The following matrix elements values were returned 
by the fits, performed either within the fragmentation approximation 00 or evaluating the 
full leading order matrix elements 0, |T3|j: 
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m is the charm mass, usually taken equal to 1.5 GeV. Figure [L] shows the CDF data from the 
Tevatron and the curves which fit them with these parameters. 


The uncertainties on these fits are certainly not smaller than a factor of two, due to the 
many systematics entering their determination: parton distribution functions (responsible for 
the difference between [|14|] and ||0), charm quark mass, a s value, higher order QCD corrections, 
etc. They could, however, even be larger. Indications in this direction come from a fit [HJ which 
makes use of PYTHIA for simulating the effect of initial state radiations from the partons before 
they collide to produce the QQ pair: this changes the slope of the theoretical predictions and 
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Figure 1: Fits to production of direct J/tJj (left) [jl2| and x (right) jTlj] at CDF. The ,J/'ip plot also 
shows by how much the Color Singlet Model underestimates the data. Also shown (below) is the J/'ip 
polarization pattern predicted jl3|] by the factorization approach. 


hence the result of the fit: 


(08 /V, ( 3 Si)) ~ 3 x 1(T 3 GeV 3 [15] 

(OfCSo)) + -^(Of( 3 P 0 )) ~ 1.2 x 1CT 2 GeV 3 
m z 
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The effect can be seen to be large, the results being significantly smaller. While such a reduction 
would actually be welcome in the light of other data which will be presented further on, one 
should however for the time being only take this as an indication of the size of the uncertainty 
which may still lay hidden in their determination. 

Other than the size of the production cross section, the polarization of the quarkonia is an 
observable with great discriminating power for the various approaches to quarkonia production. 
It can be measured by analyzing the angular distribution of the quarkonium decay products 
(muons) in its rest frame, and parametrizing it as 


da{^ —> fi + fi ) 
d cos 6 
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Figure 2: Total cross section (left) and inelasticity distribution (right) in photoproduction at HERA 
p~ 8 [] . Color octet curves with parameters fitted to Tevatron data in JT 2 [] . 


J/0’s produced at the Tevatron at large px are predicted to be almost fully transversely po¬ 
larized, i.e. a(J/ 0 ) — 1 JT 6 j, as a result of the production via gluon fragmentation into 3 S[ 8 ^ 


states being largely dominant. At smaller px : on the other hand, non-fragmentation channels 
involving l S^ and 3 Pj 8 ' become important: the J/ 0 ’s are then predicted to be produced es¬ 
sentially unpolarized in the low transverse momentum region, around px — 5 GeV 0- The 
observation of such a polarization pattern, shown in figure [T], would provide great support for 
the factorization approach to quarkonia production. 


3 The HERA data in 7 p collisions 

Within 7 p collisions ,//?//s can be produced in leading order at non-zero px via the color singlet 
channel yp —> 3 S^g —> ,J/yjX (||, hrst reference). Next-to-leading order QCD corrections to 
this channel have been recently computed JlTj, and the results have been found in fairly good 
agreement with the experimental results from the ZEUS and HI experiments: the absolute 
normalization of the total cross section agrees within the theoretical uncertainties, and the 
shape of the inelasticity distribution of the J /0 (usually denoted by z, with z = Ej/^/E 1 in the 
proton rest frame) is well described by the calculation. 

Color octet contributions to J/7 photoproduction have been investigated in leading order 
||T8| , p~9[1 . In the non-zero px region the five yp —> ( 1 5'o®, 3 <Si^, 3 Pj^)g channels contribute. 
Figure |2] shows the results for the total cross section and the inelasticity distribution with the 
inclusion of these channels: matrix elements of the order of the ones fitted to the Tevatron data 
without using PYTHIA have been employed in these plots, taking (0^’( 3 St)) = (0 ^ p ( 1 Sq)} = 
(0^( 3 Pj))/( 2J + l)m 2 = 1 x 1(T 2 GeV 3 . 

One can see from the plots how the data do not need any octet contributions: the color 
singlet channel by itself can describe them well. More than this, the octet terms evaluated 
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with the Tevatron parameters look at variance with the data, suggesting a non-universality of 
these NRQCD matrix elements. We have however seen how indications exist [HJ that the fits 
to the Tevatron data may be a significant overestimate: if we reduce the value of the matrix 
elements used in the photoproduction calculation by a factor of three, to bring them in line 
with the smaller Tevatron fits returned by using PYTHIA, the discrepancy in the inelasticity 
distribution is greately reduced. Further unaccounted for contributions, like higher orders near 
the phase space end point, could easily provide large corrections and bring the prediction in 
agreement with the data. See refs. JJJ. ^U] for a more detailed discussion about this point. 

An analysis of photoproduction data within the factorization approach has also been at¬ 
tempted in the elastic region, by fitting experimental results with the leading order prediction 
given by the octet channels 7 g —> ^Sq^, 3 P 0 , 3 P 2 ^■ This gives the result [|2T] 


(Offs 0 )) + -^(O^CPo)) ~ 2 x 10 - 2 GeV 3 [ 21 ] 
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This looks smaller than some of the Tevatron fits but, as we shall see in the next Section, in 
line with fits to fixed target data. This result should however be taken with great care, due to 
the many subtleties surrounding elastic vector meson production and to the large theoretical 
uncertainties. 


4 Quarkonia in fixed target experiments 


Fixed target experiments were of course the first to produce data on quarkonia production. 
Large normalization discrepancies between, say, J/0 production data and theoretical predictions 
based on the Color Singlet Model had been observed but, lacking a detailed understanding of 
the problem, usually dealt with by scaling the theoretical curves by very large K-faetors, of order 
ten or more. Nowadays, the FA offers the possibility to give a theoretically sound interpretation 
of these data, possibly in line with the one which seems to successfully describe the Tevatron 
data. Color octet contributions to total J/ip and 0' production in a fixed target set-up have 
been evaluated, and compared to data obtained in pN and 7 tN collisions. Fitting pN data with 
leading order short distance cross sections (a next-to-leading order calculation is in preparation 
P2||) the following values for the matrix elements have been found ||23|| : 


(OfCSo)) + -L(of( 3 P 0 )) - 3 X 10 - 2 GeV 3 
rrr 

(Oi'CS 0 )) + — 2 (Of(*P 0 )) ~ 5.2 x 10 ~ 3 GeV 3 
rn 


[23] 

23 


(13) 

(14) 


Figure shows how such values for the color octet matrix elements allows for a good description 
on the data in a wide beam energy range, whereas the singlet contribution alone clearly under¬ 
estimates them. A similar analysis has also been performed in |24|]. It is to be noted that the 
fitted values are not fully consistent with the Tevatron ones, looking at least a factor of three 
smaller (notice that the linear combination of the two matrix elements is not exactly the same): 
this can be taken as a further indication that the Tevatron fits might be an overestimate). 

These parameters fitted in pN collisions have been tested by predicting J/0 and 0' produc¬ 
tion in 7 tN collisions: they provide a fairly good description of the data, though they seem to 
underestimate them by about a factor of two 
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Figure 3: J/fy> (left) and ?// (right) production in pN fixed target experiments, compared with 

theoretical fits with the color octet contributions [p3[. 


A few more observables can offer a good handle on the validity of a given quarkonium 
production framework, namely y’s production, relative fraction of t/ds and y’s, and polarization 
of the produced quarkonium. Let us consider them in turn. 


The ratio cr(y 1 )/cr(y 2 ) has been measured in many experiments, both in pN and in nN 
collisions. An average value for nN is about 0.6, while two pN experiments give smaller values 
(0.34 ± 0.16 for E771, 0.24 ± 0.3 for E673) but with large errors. A new preliminary result, 
0.45 ± 0.2, has been presented by E771 at this Conference p5 |. 


These experimental results have to be compared with a theoretical prediction, within the 
FA, of about 0.08 |^3|. The reason for such a small prediction is that at leading order the 
gg —+ 3 P[ 1 ' > —> Xi process vanishes in the FA, and the color octet channel qq —> 3 S^ —► \J only 
contributes a small fraction. One should however consider that more color octet channels can 
contribute to both Xi and Xi production. These processes go through the octet states 
3 Pj 8 ^ and 3 D^\ and an accurate assessment of their relevance is prevented by our ignorance of 
the values of the NRQCD matrix elements weighing their transition to the observable x states. 
A very crude estimate of their effect |23], Q returns a value around 0.3, in good agreement with 
the pN data but smaller that the nN ones 1 ]. 


One more interesting observable is the fraction of J/'ip coming from y’s decays, a(J/0 <— 
x)/<j(J/0), which is of course also a part of the ratio of J/0 and y’s cross sections. The 
experimental data, both in fixed target experiments in pN and nN collisions and also in pp 
collisions at the Tevatron, gather around a central value of 0.3-0.4. The FA, making use of fits 
to Tevatron data for the %’s matrix elements and to fixed target data for J/ip ones predicts a 
value of about 0.3 |23[, thus in good agreement with the experiment. 


What’s special about this observable is that its value is strikingly different in qp collisions. 

x Of course this disagreement between the two sets of data, if found to persist and become more significant, 
would by itself be pretty puzzling. 
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This is because the leading order reaction 7 g —► xj9 is forbidden due to charge conjugation 
invariance. Indeed, no experimental data exist for cr(J/0 <— x)/cr(J/0) bi photoproduction, due 
to the vanishingly small y’s yield, but only an upper limit by NA14 |2fJ: <- < 

0.08. The reason why this is important is that the Color Evaporation Model, in that it is 
not concerned with the details of the particles initiating the reaction, would predict the same 
ratio for hadro- and photoproduction. Hence, one could say that the CEM is ruled out by 
this result, though more independent confirmations of the photoproduction experimental result 
would actually be welcome. 

One further key observable is the polarization of the quarkonium, measured according to 
eq. ( 0 )- One finds, experimentally, 


a( J/ip) = 0.02 ± 0.14 (15) 

a(V0 = 0.028 ± 0.004 (16) 


That is, the quarkonia are found to be produced essentially unpolarized. 

This is in contrast with theoretical calculations within the FA, which predict instead a 
sizeable degree of polarization, returning 0.31 < a(J/yj) < 0.62 and 0.15 < a('0 / ) < 0.44 


the large band taking into account the very approximate knowledge we have of the NRQCD 
matrix elements’ values. 


Such a discrepancy is certainly disturbing, and if conhrmed would be a serious problem 
for the factorization approach. At the present stage of our understanding we must however 
be aware that many unaccounted for contributions may still play an important role here. For 
instance, a proper inclusion of effects which would lead to off-shell rather than on-shell colliding 
gluons could significantly change the picture presented above. I therefore think it would not 
be wise to bury the factorization approach at this stage and because of this discrepancy. 


5 Conclusions 

With this brief survey of a few of the phenomenological consequences of the factorization 
approach to quarkonia production (which I should urge you not to call “Color Octet Model”!) 
we have seen how it looks able to describe in a satisfactory way experimental results previously 
at great variance with the Color Singlet Model, both in fixed target experiments and at the 
Tevatron. 

Problems appear in photoproduction at HERA, where the color singlet contribution alone 
appears on the other hand to well describe the data. But the uncertainties are still large and 
can accomodate for the discrepancy. 

Polarization data are also troublesome for the factorization approach (as they are for the 
CSM, while the CEM cannot give a prediction at all). But I feel, once again, the uncertainties 
to be still too large for an assessment of the validity of this approach based on these data. 

Surely enough, the factorization approach appears superior to both the CSM (of which it is 
an extension) and the CEM. One could say the CSM approximation still appears to work fairly 
well when few gluons are involved (in photoproduction), whereas the CEM can work when very 
many gluons are around (in hadron-hadron collisions). But neither of them can even attempt 
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to describe the whole yield of data: for instance, the CSM can be ruled out by the Tevatron 
data alone, and the CEM could possibly be ruled out by \j photoproduction. 

It is quite a widespread belief (though by no means universal!) that the factorization 
approach is the right theory for heavy quarkonia production and decay. As a matter of fact, 
the NRQCD lagrangian on which it is based is nothing but a limit of the QCD lagrangian itself, 
and not an ad hoc model. 

If anything, a problem of the FA is not of being inadequate, but rather of being perhaps even 
too general. Many different NRQCD matrix elements enter the phenomenological predictions 
(because the corresponding operators enter the nonrelativistic limit of the QCD lagrangian), 
and it is difficult to produce accurate numerical results with so many unknown parameters. 
These matrix elements are on the other hand rigorously defined and in principle calculable by 
lattice QCD, so it is possible they will be more precisely determined in the future either this 
way or by global analyses of experimental data. 

Given the reasonable correctness of the underlying lagrangian, discrepancies of the theoret¬ 
ical predictions with experimental data may still be originated by the approximations included 
in the calculations. 

Higher twist corrections to the factorization approximation can be large, especially for 
charmonium (one should never forget that A jm c ~ 0.3, not really a negligibly small number). 
Higher order corrections, both in the strong coupling ( a s {m c ) ~ 0.3) and in the velocity of 
the heavy quarks (again, v 2 ~ 0.3 for charm), can also be large. The value of v 2 (and higher 
powers) determines - via scaling rules [^7]] - which matrix elements are dominant, and also how 
accurately - via spin symmetry - different matrix elements can be equated to each other. These 
approximations are widely used in the phenomenological predictions, to truncate the series and 
to decrease the number of independent parameters. How reliable they are is therefore extremely 
important for the accuracy of the theoretical result: the fairly large value of v 2 for charmonium 
systems may help explaining discrepancies between nowadays predictions and experimental 
data. 

The situation should be much better for bottomonium systems, for which all the expansion 
parameters I mentioned take significantly smaller values, around 0.1. All the approximations 
should therefore be much better justified, and one should expect a better agreement between 
theory and data. A detailed study of such systems, both theoretically and experimentally, will 
therefore greatly help finally confronting the factorization approach to quarkonia production 
and decay with the real world. 
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